Theory of many-fermion systems II: The case of Coulomb interactions 
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In a recent paper a general field-theoretical description of many-fermion systems with short-ranged 
interactions has been developed. Here we extend this theory to the case of disordered electrons 
interacting via a Coulomb potential. A detailed discussion is given of the Ward identity that 
controls the soft modes in the system, and the generalized nonlinear a model for the Coulombic case 
is derived and discussed. 



I. INTRODUCTION 

The theory of many-fermion systems is central to our 
understanding of condensed matter, as well as of nu- 
clei and certain astrophysical systems. Of particular in- 
terest in a condensed matter context are the universal 
properties of many-electron systems, i.e. phenomena at 
long wavelength and small frequency scales that are inde- 
pendent of the material's detailed microscopic structure. 
Historically, there have been two important avenues to 
this priiblem: Landau's phenomenological Fermi-liquid 
theory,El and the microscopic many-bgdy perturbation 
theory or Feynman diagram approach.l3 The latter soon 
was generalized to include the scattering of electrons by 
static impurities ,13 and it has led to many important in- 
sights concerning the nature of interacting disordered 
electron systems. For instance, it was used to show that 
the combined effects of disorder and interactions lead to 
the nonanalyticities in the frequency and wavenumber 
dependence of both thermodynamic and transport prop- 
erties of disordered metakjj-that have become known as 
'weak localization effects'.ffl 

There are, however, many interesting phenomena for 
which an approach based on many-body perturbation 
theory is not feasible. An example is the Anderson- 
Mott metal-insulator transition that the electrons un- 
dergo with increasing disorder strength.l3 This quantum 
phase transition is best studied by means of effective field 
theories and the use of the renormalizatipn group, an 
approach that was pioneered by Wegner.Q Also, many- 
body perturbation theory is ultimately unsatisfactory as 
a tool even in the metallic phase, despite its impressive 
successes. One major problem is that, within perturba- 
tion theory, it is not clear whether the weak-localization 
effects are actually the leading nonanalyticities. Further- 
more, universal phenomena like the weak-localization ef- 
fects arise from the presence of soft modes in the problem. 



which many-body perturbation theory is not well suited 
to deal with. The softness or masslessness of these modes 
is the result of symmetries that are often not explicit 
in the usual perturbative formalism, and therefore soft 
modes arise apparently accidentally as a result of com- 
plicated cancellations, rather than being manifestly built 
into the formalism. Finally, it is unsatisfactory to have 
entirely different approaches to the universal properties 
of the various phases on the one hand, and to those of 
the transitions between these phases on the other hand. 

The concepts for overcoming these problems exist; they 
consist of a systematic application of renormalization 
group ideas. Historically, the renormalization group has 
been mainly associated with critical phenomena at con- 
tinuous phase transitions, for which the importance of 
symmetries and soft modes has been universally appre- 
ciated. While it has always been clear in principle, and 
occasionally has been emphasized,El that the renormal- 
ization group, rather than just being a tool for studying 
critical phenomena, allows for a unified description of 
both phases and phase transitions, these ideas have not 
been widely appreciated. Only very recently has there 
been the beginning of a paradigm shift in this respect. 
For instance, starting with Shankar's work,la there has 
been much activity recently on the derivation of clean 
Fermi-liquid-theory as a stable renormalization group 
fixed point.E3 These methods, however, have proven very 
hard to generalize to the case of quenched disorder. 

In a previous paper, to be referred to as I,lliJ two of us 
have developed an effective field theory for many-electron 
systems that is particularly suitable for dealing with dis- 
ordered systems. This theory, which is formulated in 
terms of classical matrix fields, allows for a systematic 
separation of soft and massive modes, and the latter can 
be integrated out in a simple approximation to yield an 
effective theory for the soft modes. The theory also al- 
lows for the clean limit to be taken, but the soft-mode 
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structure in that case turns out to be more complex, 
which gives the effective theory fewer advantages over 
traditional approaches than is the case in the presence 
of disorder. This theory has been used, inter alia, to 
prove that the well-known weak-localization effects are 
indeed the leading nonanalyties, to provide a technicaliy, 
satisfactory derivation of Finkel'stein's generalizationta 
of Wegner's nonlinear a modelJ3 and to derive a previ- 
ously unknown nonanjalvlicity in the spin susceptibility 
of clean Fermi liquidsJljO 

In I, the effective classical field theory was developed 
for fermions that interact via a short-ranged potential, 
assuming that the underlying Coulomb potential had 
been screened at the level of the basic fermionic theory. 
In the present paper, we show how to handle a long- 
ranged interaction entirely within the framework of the 
classical matrix field theory, focussing on the disordered 
case. We give a detailed discussion of the Ward iden- 
tity that controls the structure of the soft modes, apd we 
derive Finkel'stein's generalized nonlinear a modelij for 
the Coulombic case. 

The outline of our paper is as follows. In Sec. ^ we 
recall the matrix field theory of I, and slightly generalize 
it to allow for a Coulomb interaction. We show that the 
saddle point considered in I remains valid, and that an 
expansion about the saddle point to Gaussian order pro- 
duces RPA-type screening. In Sec. Ill we perform a sym- 
metry analysis, and derive and discuss the Ward identity 
that controls the soft modes in the system. In Sec. IV 
we show that integrating out the massive modes in the 
simplest approximation that respects the Ward identity 
leads to Finkel'stein's model. In Sec. we conclude by 
discussing our results, and give in particular a discussion 
of the accuracy of the a model. 



II. MATRIX FIELD THEORY 
A. Grassmannian field theory 

Our model and basic theoretical setup is the same as 
in I. We will therefore restrict ourselves to discussing the 
changes that are necessary to accomodate a Coulomb in- 
teraction. 

The action is given by 



S ■ 



dx ^ ■0(7 (x) dr ijja- (x) 



So + Sd 



(2.1a) 



Here the ip and ip are Grassmann valued fields, and 
we use a (d -|- l)-vector notation, with x — (x, t), and 

J dx = Jy dx Jq dr. x denotes position, t imaginary 
time, V is the system volume, P — 1/T is the inverse tem- 
perature, a is the spin label, and we use units such that 
h — kB = = 1- So (together with the time derivative 



term) and S'dis describe free fermions, and their inter- 
action with a static random potential, respectively, and 
have been defined in I. The random potential we assume 
to be Gaussian distributed, and we employ the replica 
trick to handle it. Sint describes a Coulomb interaction, 
which in the replicated theory (which we denote by a 
tilde) takes the form 

^int = -^int = ^\Y j '^^l X! "(^^ ^ ^2) 



X5(ri - T2) < [X,] < {X2) (^2) V'", (^1) 



Here a is the replica index, and 
«(x) = l/|x| , 



(2.1b) 



(2.1c) 



is the Coulomb potential. For the dimensionalities of 
interest, d = 2, 3, its Fourier transform is 



i'(q) = (1 - 6c^fi) 



2'^-V 

\r^\d-l 



(2.1d) 



where the factor 1 — (5q^o represents a uniform positive 
background charge that ensures charge neutrality. 

We now introduce a momentum cutoff A, and rewrite 
the interacting part of the action as 



qa _ c"(l) ( c"(2) _L c"(3) 
•^int ~ *^int *-'int "-"int 



(2.2a) 



where, 



S" 



't '^ = - 2 ^ ^ ^ '"'('^^ + '^^ 

(71(72 k^p q 



(2.2b) 



S; 



a (2) 



T 



EEE'^(p-k)0(ip-ki-A) 



o"i(72 k.p q 

xC,(fc)V;?Jp + g)C.(fc + 9X(p) , (2.2c) 



^in/'^--| E EEVk + p)e(|p + k|-A) 

xC, {-k) C {k + q) V4 (-P + 9) C b) • 

(2. 2d) 

As in real space, we use a (d -I- l)-vector notation with 
k = (k, cjn), where ujn — 2'KT{n + 1/2) is a fermionic 
Matsubara frequency. The prime on the g-summation 
indicates that only momenta up to the momentum cut- 
off A are integrated over. While the long- wavelength, 
small-frequency phenomena we are interested in do in 
general not depend on A, the choice of this cutoff will be 
important for the range of validity of the final effec tive 
theory. We will come back to this point in Sees. HI and 
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^ below. Equations ( |2.2| ) represent the same phase space 
decomposition as in Eqs. (2.8) of I, except th at we hav e 
expHcitly inserted the step functions in Eqs. ( 2. 2c , 2. 2d ), 
because the s mall- w avenumber part of v is already con- 
tained in Eq. (2.2b). In I the overcounting that resulted 
from not having the step functions explicitly present was 
of no consequence, since v in that case was not singular in 
the small- wavenumber limit. For a Coulomb interaction, 
more care must be taken. 
We next introduce spinors 



and their Fourier transforms 



(2.3a) 



^"(fc)^V"(k)= (^^ijjj) , (2.3b) 

as well as their adjoints, ■0"(fc), and a scalar product in 
spinor space, (tpTip) = "0 ' V'l where the dot denotes the 
matrix product. Then we can write the interaction term 
as 



qa _ («) I (*) ( 5" (3) 
"^int — '^int *-'int '-'int ' 



(2.4a) 



with, 



k,p q 

x(^"(p + g),so^"(p)) 



(2.4b) 



that were defined in I. The result of this procedure, which 
was explained in Appendix A of I, is 



k,p q 

x{^'-{p + q),s^r{p)) , (2.6a) 



I? 



a(t) 



"2 



■r^*^EE E(v^"(^)'^^^"(^ + '?)) 

k,p q i—1 

x(V'"(p + (z),5.^"(p)) . (2.6b) 



5 a (3) 
'-'int 



T _ 
2 



■r^'^ E EE c.(fcX(-^+9) 

<7iyia2 k,p q 

xC.(P + 9)Ci(-p) ■ (2.6c) 



Here 



r(^)(q)=z;(q)-r(*) 



(2.6d) 



and r(*) and F^^' are numbers t hat r esult from integrating 
over the wavevectors in Eqs. ( |2.5| ) as explained in Ap- 
pendix A of I. Notice that as a result of this procedure, 
r^*) and r'^'^^ for clean electrons depend logarithmically 
on the cutoff A, and diverge as A ^ 0. For the disordered 
case, the logarithmic singularity is protected both by A 
and by the disorder. This singularity, which is a conse- 
quence of the Coulomb interaction, is the reason for our 
modification of the procedure employed in L 



srJ' = E E 'rg{q)J2{r{k),s.r{k + ?)) 

k,p q i—1 

x(V'"(p-t-g),5.^"(p)) . (2.4c) 

Here Sj — iaj, with (7i^2,3 the Pauli matrices, and sq = ctq 
is the 2x2 identity matrix. We have also defined the 
singlet (s) and triplet (t) interaction amplitudes 

r£('?) = ^«(p-k)e(|p-k|-A) , (2.5a) 

and 

rS(<z)-^'(q)-ri*],(g) . (2.5b) 

In addition we define the Cooper channel or 2k p- 
scattering amplitude 

^g(q) = ^;(k + p)e(|p + k|-A) . (2.5c) 

These expressions are the same as the corresponding ones 
in I, with the exception of the momentum restrictions 
in the effective interaction potentials r^''^''^^ discussed 
above. Next we project the modes in ^int onto density 
modes in the particle-hole and particle-particle channels 



B. Composite variables: Matrix field theory 

After completing the phase space decomposition and 
projecting onto densities we are in a position to refor- 
mulate the theory in terms of composite variables. This 
proceeds in exact analogy to I, and we therefore only 
quote the result. The partition function Z for the repli- 
cated theory is written as an integral with respect to two 
matrix fields, Q and A, 

Z = J D[Q] D[A] , (2.7) 

with an effective action 

A[Q, A] - AdrsiQ] + AntiQ] + i Tr In (^G^^ - zA) 

+ J dxtY (A(x)g(x)) . (2.8) 

The matrix field Q corresponds to expressions that are 
bilinear in the fermionic fields ip and ^. Correspondingly, 
Q carries two frequency indices n and to, and two replica 
indices a and /3. Each matrix element Q^m is an element 
of Q X Q, with Q the quaternion field. A is an auxiliary 



3 



field whose technical role is to constrain the products of 
fermionic fields to the Q. It is convenient to expand the 
in a spin-quaternion basis, 



- nm 



Q 



af3 
nm 



rA—O 



a/3 



(2.9) 



and analogously for A. Here tq = sq = I2 with I2 the 
2x2 unit matrix, and Tj — — Sj = —ioj (j = 1,2,3), 
where the Oj are the Pauli matrices. The properties of 
Q and A have been derived in I, and for c omp leteness 
we list them again in Appendix In Eq. (2^) and in 
what follows, Tr is a trace over all degrees of freedom, 
including an integral over x, while tr is a trace over all 
discrete degrees of freedom that are not shown explicitly. 

Gq-i = -9^ + VV2m + ^ , (2.10) 

is the inverse of the free electron Green operator, a nd i t 
is clear from the structure of the Tr In-term in Eq. ( |2.§| ) 
that the physical interpretation of the field A is that of a 
self-energy. The contributions Adis and ^int to the action 
read 



(2)r 



-1 



2'kNfti 



ctx. 



(trQ(x; 



[Q] = 



dx tr(Q(x))^ 



(2.11a) 



(2.11b) 



(2.11c) 



with Tioi the single-particle relaxation time and ti a re- 
lated scattering time defined in I, and 



■Aint [Q] — -^int + •^int + 



(2.12a) 



A 



(s) 



r— 0,3 ni ,712 ,m a 

X [tr {{Tr^So)QZn,+mi^))] 

^/dxE(-ir E EE 

r— 0,3 ni,n2,m a i—1 

X [tr ((r.®S.)Qnr.n.+mW)] 

X [tr ((r,®s.)Q^Jr+™,„,(x))] , 



(2.12c) 



^-E E E 



r—1.2 ni ,712 ,771 a 
X [tr {{rr^So)QZ-rn+rni^))] 

X [tr ((r,®so)Q^"„,,„,+,„(x))] . (2.12d) 

Here Np is the density of states at the Fer mi leve l in 
saddle-point approximation, as defined in Eq. ( ^.13c ) be- 
low. We have written the action in real space, but one 
should remember that all of the fields are restricted to 
Fourier components with wavenumbers |k| < A. 



C. Saddle-point solution, Gaussian approximation, 
and physical correlation functions 

It is easy to see that the Fermi-liquid saddle-point dis- 
cussed in I remains a valid saddle-point in the presence 
of a long-range interaction, with F^'*^ in I replaced by 
F(")(q = 0). The single-particle Gr een function in saddle- 
point approximation is therefore given by the same ex- 
pression as in I, viz. 



G'sp(p, uj„) = [iuj„ - p^/2to + h - Sn] 



(2.13a) 



with fi the chemical potential, and the self-energy S the 
solution of the equation 



— 



V ^ 
p 



lijJr, 



p2/2TO 



(2.13b) 



This is the standard Hartree-Fock result, with the disor- 
der treated in self-consistent Born approximation. 



-2 1 



A^F = — ^ VlniG,p(p, 



iO) 



(2.13c) 



is the density of states at the Fermi level in saddle-point 
approximation, which is used for normalization purposes 
throughout. 

Similarly, in an expansion about the saddle-point to 
Gaussian order, all of the expressions derived in I remain 
valid if we substitute r*^*)(p) for F^*) in all propagators 
taken at wavevector p. In particular we obtain for the 
low-frequency, long-wavelength limit of the density sus- 
ceptibility, X, in the disordered case 



X(p,f^„) = Xst(p) ■ 



d(p)p2 



|a.|+d(p)p2 ■ 
Here r2„ is a bosonic Matsubara frequency. 

-Np 

xst(p) = i^^^r(.)(p) ' 

is the static density susceptibility, and 
d(p) = i?[l + 7VFF(^)(p)] , 



(2.14a) 



(2.14b) 



(2.14c) 



with D the Boltzmann diffusion coefficient. In the clean 
limit we obtain the usual RPA expression, 



x(p,r2„) = 



Xo(P,f^n) 



i + r(-Hp)xo(p,f^„) 



(2.15) 



with xo the Lindhard function. We see that in Gaus- 
sian approximation, the field theory describes screening, 
and the existence of plasmons, in RPA and its disordered 
generalization, respectively. 
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III. SYMMETRY ANALYSIS 



From Sec. II C it follows, in conjunction with I, that 
at the level of the Gaussian approximation the Qnm with 
nm < are soft modes. In this section, we perform a 
symmetry analysis in order to show that they are indeed 
the exact soft modes of the theory. This will allow us 
to explicitly separate the soft modes from the massive 
ones, and to formulate an effective field theory for which 
the soft modes remain manifestly soft to all orders in 
perturbation theory. 



A. Basic transformation properties, and Ward 
identity 

Let us start with a symmetry analysis of our field the- 
ory that is a slight generalization of the procedure fol- 
lowed in I, which in turn was a generalization of the work 
on noninteracting eleciipns by Schafer and Wegner, and 
Pruisken and Schafer .Efl We consider an infinitesimal si- 
multaneous rotation in frequency and replica space space 
given by. 



nm 



(3.1) 



which transforms the Q-matrices according to Q — > 
TQT~\ with T = CfC^, where C = in ® sa- {ni > 
0, n2 < 0) and {ai,a2) are fixed pairs of frequency and 
replica indices that characterize the transformation. Un- 
der such an infinitesimal rotation, the Q-matrices trans- 
form like 



jaf3 
i nm 



a/3 
nm 



with 



)Q/3 ^ 

^ nm 



a2p 

712 m 



^/3ai ^mn 



1 Vnn2 



-(1^2)) + o{e^) 



(3.2a) 



(3.2b) 



Here we have shown only the frequency and replica in- 
dices, since all other degrees of freedom are unaffected 
by the transformation. The A-matrices transform ac- 
cordingly. The symbol (1 ^ 2) indicates the same terms 
as written previously, but with the indices 1 and 2 inter- 
changed. 



Of the terms in the action, Eq. (2^), only Tr ln(GQ ^ — 
iA) and Aint are not invariant under the above transfor- 
mation. Their transformation behaviors are easily deter- 
mined by an explicit calculation. We find 

Tr In (Go^i - iA^ ^ Tr In ^ - iA^ + 6 Tr (G Situ) , 

(3.3a) 



withG=(Go^-iA)-\and 

li^^^)nm ~ ^iO^rO {^aai^nrii^ fja2^mn2 + Saa2^nn2 

X^;3ai^rrmi) if^iii-na : (3.3b) 



and 



As) ^As) ,XAis) 



(3.4a) 



with 



<5A1:]=32 /'dxdyr(^)(x-y)^ 

r=0,3 

x^^E[°C:w?C:U(n„-«.)(y) 



-(1 



(3.4b) 



For our purposes, we concentrate on a discussion of the 
particlc-holc spin-singlet interaction; the discussion of 
the remaining interaction channels proceeds analogously, 
but is less interesting since in these channels the interac- 
tion is short-ranged. Proceeding analogously to I, we ob- 
tain from the above transformation properties the Ward 
identity 



W,nt + 8nn,-n2 J dy {''oQM oQsii^)) 
= <5l3<S24((|]Qn(x))-(E!Q22(x))) , 



(3.5a) 



where, 



W^nt = -32 f dydz r(^)(y - z) ^ T 

r=0,3 

(^) 



r— 0,3 nanb 



-(1^2) 



(3.5b) 



To simplify the notation we have adopted the convention 
1 = {ni,ai). For the special case of a short-ranged in- 
teraction, r(*)(x — y) = r^'''' (5(x — y), we recover Eqs. 
(3.14) of I. 



B. Solution of the Ward identity 

We now solve the Ward identity by employing a 
method that is more transparent than the one used in I, 
and more suitable for a generalization to the long-range 
case. 

We first split up the Q-matrices into their averages and 

fluctuations: Qnm = (Qnm) + (A(5)„m = 6nm{Qnn) + 

(A(3)„,„. Since ni 712, and ^ 714, we have (Qmn-,) = 
{Q nsrii) — 0. Furthermore, if we put Ua — nh in the 3- 
point functions, then the expression in square brackets 
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in Eq. ( |3.5b| ) vanishes due to Eq. ( A3a ), so effectively we 
also have (Qn^nt) = 0. Equation (|3.5a ) then takes the 
form, 



Wint + 8r2„^_„2Cl2,34 — Sl3 S24 Ni: 



(3.6a) 



where 



Ci2,34 = y'dy ^0(AQ)i2(y)0(AQ)34(x)^ , (3.6b) 



and 



(3.6c) 



Wint can be written as the sum of two terms, Wir. 

W^inV+l^iL't\with 

.(1) 
' int 



W>„V = -32(5„,„,r(^)(k ^ 0)iVi2 



(3.6d) 



and 



r=0,3 q riaUi, 

' °(AQ)::";,_(„„_„,)(-q) "AmzTM) 

xS](AQ)34(x))-(1^2) , (3.6e) 

(2) 

Here we have chosen a mixed representation for W-^^^ , 
with some AQ in real space, and some in Fourier space, 
in order to make the cutoff on the q- integration explicit. 

To proceed, let us ignore W^\^^ for th e tim e being. Its 
effect will be analyzed later. Equation (3.6a) then turns 
into a closed integral equation for the homogeneous corre- 
lation function C. At this point we note that our global 
symmetry transformation, Eq. (3.1), produces a Ward 
identity for homogeneous correlation functions. For a 
short-ranged interaction, this is sufficient to capture the 
important structural restrictions imposed on the theory 
by the symmetry of the action. However, in the long- 
range case the homogeneous limit is singular, and we have 
to be more careful. For instance, it is obvious that a lo- 
cal symmetry transformation would generate a wavenum- 
ber dependent P'^'*^ in VK/jj^^^ and in the long-range case 
r(«)(k ^ 0) 7^ r(^)(k 0). This is the reason why we 
have written r('^^(k 0) in Eq. (3.6d). Furthermore, it 
is known from perturbation theory that the dispersion of 
the soft modes that_are controlled by the Ward identity 
is diffusive, ~ k^.E£l The equation for the wavenumber 
dependent C then takes the form 



+ Dk"^) Ci2,34(k) = (5i3 ^24 NuCk) 
+ 32r(^)(k) A^i2(k) 

^l-2,a-b C'a6,34(k) , 

ab 



(3.7) 



with D the (exact) diffusion constant, and 7Vi2(k) the 
wavevector dependent generalization of A'^12 as defined in 
Eq. (3.6c). It is useful to define 



-fri(k) — T ^ ^ Sn.ni—n2 -^nin2(k) 



(3.8a) 



By performing appropriate summations over the Ward 
identity, wc find that 



p„(k) = i (r!„ + Dk2) xsc(k,rj„) 



(3.8b) 



with Xsc the screened density susceptibility. The latter 
is defined as 



Xsc(k, Qn) = 



i + r(-)(k)x(k,f]„) ' 



(3.8c) 



where 



X(k,r!„) =32r^^„,„,_„,Ci2,34(k) , (3.8d) 



1234 

is the full density susceptibility. Notice that 
this relation between P„ and Xsc is exact, since 

X/12 '^ri,"l-n2^12,34 ~ ^' 

The integral equation, Eq. (3.7), can be solved in terms 
of P„ or Xsc by means of the same methods that were 
employed to discuss the Gaussian Q-field theory in I. We 
find 



Cl2,34(k) = SisS24'Di2{k) 

16 



+^l-2,3-4<5aia3<5aia2 2TT*^'*-'(k) 

xl?i2(k)2?(:^(k) 



(3.9a) 



where 



I?i2(k) =2iVi2(k) [au-n^ +^k2] ' , (3.9b) 



and 



v[l\k)=Vi2{k) i-r(^)(k)xsc(k,r!„,-„j 



(3.9c) 



To appreciate the difference between this structure and 
the analogous one in the short-range case, it is instruc- 
tive to consider the limit of small wavenumbers and 
frequencies with Dk^ r2„j_„2. In this limit we 
have 2iVi2(k) TTN{eF) with N{eF) the exact den- 
sity of states at the Fermi level, and Xsc(k, 57,n_„2) ~ 
k^/r2„j_„2- In rf = 3, we obtain 



16Ci2,34(k^0) = 512^: 



34 



TrNiep) 

—712 



o nk2 
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We see that the long-ranged Coulomb interaction causes 
the part of C that is non-diagonal in frequency to diverge 
like T/r2k^, rather than like T/fi^ in the short-range 
case. This structure is responsible for the well-known log- 
squared terms in the density of states and in the sound 
attenuation coefficient that appear in perturbative cal- 
culations of weak-localization effects in 2-d systems with 
Coulomb interactions.Q 

We now consider the remaining contribution to Wint, 

(2) 

W^ijjj . Without the restriction on the wavenumber inte- 



gral in Eq. (|3.6q ), one could relabel wavevectors to show 
that VFjnt contains terms that have the same structure 
as W^int^. However, due to the cutoff contained in the def- 
inition of our Q-field theory this is not the case, and we 
resort to perturbation theory to analyze the structure of 



(2) 



By re-expressing AQ in fermion fields, and using 



Wick's theorem to write the correlation function in Eq. 
( |3.6eD in terms of Green functions, we find 



<'t^~r(^nk)^Min(A^ 



A- 



(3.11) 



with A the cutoff introduced in Sec. [I A , and £ the elastic 



scattering mean-free path. Based on this information, we 
now choose the cutoff. There are three obvious possible 
choices: One can make A a fixed, small fraction of either 
the Fermi wavenumber kp, or the mean-free path or 
the screening wavenumber k. The first choice would not 

allow us to control the terms in VFj^^j , and is therefore 
undesirable. The second choice is possible, but means 
foregoing the option to take the clean limit, for reasons 
discussed after Eq. ( p.6d ). We therefore choose the third 



(2) 

option, which makes W^^^ of higher order in the interac- 
tion than the other terms in the Ward identity. To linear 
order in the interaction the solution given above in Eqs 
(^!9|), ( |3.10 ), is then exact. To higher order, while W'^^^ 
cannot make the modes Qnm {nm < 0) massive, it will 
change the prefactors of the diffusive singularity at small 
frequencies and wavenumbers. We conclude that an effec- 



(2) 



five theory that respects Eqs. (3.9) will have the correct 



soft-mode structure as it follows from the symmetry of 
the action. It will further exactly reproduce perturbation 
theory to first order in the interaction. In Sec. |^we will 
show that the generalized nonlinear a model of Ref. [ij 
is an effective theory with these properties. 



C. Separation of soft and massive modes 

From the previous subsection we know that the corre- 
lation functions of the Qnm with nm < remain soft, 
while those with nm > remain massive even in the 
presence of long-range interactions. Therefore the mode 
separation for this case is the same as it is for short-range 
interactions, and we can restrict ourselves to summariz- 
ing the results of 1. 



One consequence of the symmetry properties of the 
Q-matrices (see Appendix ^ is that the set of Q is iso- 
morphic to the set of antihermitian 8Nn x 8Nn matrices. 
As shown in I, from this it follows that the most general 
Q can be written as a product, 



S ps- 



(3.12) 



Here P is a matrix that is blockdiagonal in frequency 
space, 



P> 
P< 



(3.13) 



with elements P> {P"^) for nm > {nm < 0) that 
are isomorphic to the set of antihermitian ANn x ANn 
matrices, and S an element of the homogeneous space 
USp(8iVn,C)/USp(4Arn,C)xUSp(4A^n,C), i.e. the set of 
ah cosets of liSp(8A^n,C) with respect to USp(4iVn,C) x 
USp(4iVn,C)M 

This achieves the desired separation of our degrees 
of freedom into soft and massive ones. The massive 
degrees of freedom are represented by the matrix P, 
while the soft ones are represented by the transforma- 
tions S £ USp(8iVn,C)/USp(4Afn,C) x USp(4iVn,C). 

In order to formulate the field theory in terms of the 
soft and massive modes, one also needs the invariant mea- 
sure I[P], or the Jacobian of the transformation from the 
Q to the P and the S, defined by 



J D[Q] ... = j D[P]I[P] J D[S] 



(3.14) 



We will not need the measure explicitly for our purposes, 
and refer to I, where it has been constructed in terms of 
the eigenvalues of P. 



IV. EFFECTIVE FIELD THEORY FOR 
DISORDERED INTERACTING FERMIONS 

Having achieved a separation of soft and massive 
modes, we are now in a position to formulate an effec- 
tive theory for electrons with a long-range interaction 
that focuses on the soft modes. In the short-range case, 
this was done by integrating out the massive modes in 
tree approximation. This led to the nonlinear a model 
in I, and it was shown that this procedure preserves the 
structure of the Ward identity. As a result, the cr model 
contains the same Fermi liquid fixed point, as well as the 
leading corrections to scaling near it, as the underlying 
full model, and it also contains a critical fixed point that 
describes an Anderson-Mott metal-insulator transition. 

We cannot simply repeat this procedure for the present 
case of long-range interactions, since integrating out the 
massive modes in tree approximation would lead to a the- 
ory that violates the Ward identity. We therefore must 
treat the massive modes more carefully, and our aim is to 
find the simplest approximation that will still guarantee 
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the correct structure of the Ward identity, and hence lead 
to an effective theory that has the correct symmetry. 

The first steps are the same as in I: We define a trans- 
formed auxiUary field A by 



A(x) =5-i(x)A(x)5(x) , 
and a new field Q by 

g(x) = -l-5(x)(P)5-i(x) 



(4.1a) 



(4.1b) 



We then expand about the expectation values of A and 



P = {P)+AP , A = (A) + AA 



(4.2) 



As explained in I, it is sufficient to replace (P) and (A) 
by the respective saddle-point value s, and to further re- 
place (P) in the definition of Q, Eq. ( 4.1b ), by the simple 
approximation 



with 



{P)^-Np7: 



7ri2 = Si2 sgnwni 



(4.3a) 



(4.3b) 



We mention that there is no obvious small parameter that 
controls these approximations. Rather, they will be jus- 
tified a posteriori by the fact that the resulting effective 
field theory, t he n o nlinea r a model, respects the Ward 
identity, Eqs. ( |3.9[ ), (3.10), that was derived in the previ- 
ous section. This in turn shows that the approximations. 



Eqs. (4.3), are consistent with neglecting the term W-^ 



(2) 



in the Ward identity, which itself is perturbatively con- 
trolled for small interaction strengths. If the theory is 
renormalizable, this implies that the effective theory re- 
sulting from the above approximations will have the same 
structure as the full one, albeit with different coefhcients. 
We will come back to this point at the end of this section, 
and in Sec. ^ below. With these approximations, Q has 
the properties 



Q"(x) = l, QT^g, trg(x) = 0, 
and can be parametrized by 



- qiq 



(4.4a) 



(4.4b) 



where the matrix q has elements (j„,„ whose frequency 
labels are restricted ta_p > 0, m < 0. S can also be 
expressed in term of qf3 



S = 



Vl - bb^ b 

-6t VI 



(4.4c) 



1 



with 



fix) = J- (l-^/^^) 



1/2 



(4.4d) 



(4.4c) 



Unlike I, where we just dropped the fluctuations of P 
and A, here we next expand to second order in AP and 
A A. The reason for this change of procedure compared 
to the short-ranged case is that here the AP fluctua- 
tions are multiplied by a divergent Coulomb potential, 
and therefore must be retained. The part of the action 
that is quadratic in these massive fluctuations reads 



i J dxdytrG,p(x-y)AA(y)G,p(y-x) AA(x) 

-I- J dxAA(x) AP(x) 
+aI'^,[S{P)S-^ +SAPS-^] , 



(4.5) 



with A y^l \Q] the spin-singlet interaction term from Eq. 



( p. 121 ). Adis, and the invariant measure, expanded to 



(2) 

second order in AP, also contribute to Am . However, 
their net effect is to add a constant to the singular in- 
teraction potential in the A[^^ contribution to Am"^ , and 
hence they can be neglected. Likewise, there are terms 
linear in AA that couple to the soft modes S. These al- 
ways contain at least one frequency or gradient squared, 
and therefore are unimportant for the leading structure 
imposed by the Ward identity. They also turn out to be 
renormalization group irrelevant at the Fermi liquid fixed 
point that we will discuss in Sec. ^ below. We therefore 
neglect all of these terms. 

Because of the coupling between AP and S, Am still 
represents a compli cated quadratic form. To handle it, 
we expand S, Eq. (4.4c), in powers of g. To lowest or- 
der, we just have 5 = 1. It turns out that this lowest 
order approximation is sufficient to ensure the correct 
structure of the Ward identity. We have also explicitly 
checked that higher order terms in this g-expansion lead 
to corrections that are irrelevant near the Fermi liquid 
fixed point. With S — I, the massive Gaussian fluctua- 
tions are easily integrated out. Neglecting terms that are 
of first or higher order in the external frequency, the re- 

(s) 

suit is a change of the interaction term ^j^^^ to a term of 
the same structure, but with F^*' replaced by its screened 
counterpart 



F(^)(p) 



1 + NfT(^Kp) 



(4.6) 



where 



with Np from Eq. (2.13c) (here we have neglected a 
wavenumber dependence that is subleading compared to 
that of F^'')). We see that integrating out the massive 
fluctuations in the approximation we have chosen leads 
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to static screening of the Coulomb interaction. Analo- 
gous screening effects occur in the remaining interaction 
channels. However, they are uninteresting there since 
they just renormalize the numbers F*^*-* and T^'^\ 

The remaining steps in the derivation of the nonlinear 
a model are the same as in I. We thus obtain the a model 
action 



A 



NLcrM 



dx tr ( VQ(x 



-1 
2G 



+2H / rfxtr (l^Q(x) +Aint[Q] 



(4.7) 



where Q ^ Q — n, with vr from Eq. ( 4.3b ). G ~ S/ttctq 
with (To the conductivity in self-consistent Born approx- 
imation, and H — irNp/S. Aint is given by Eqs. ( [2.12 ), 



but with r(^) replaced by , Eq. (|J). This is the 
generalized nonlinear a model for disordered electrons 
with a CoulomhJnteraction, as proposed and discussed 
by Finkel'steinjlj and the above procedure represents a 
technical derivation of this model. Its properties have 
been reviewed and discussed in detail in Ref. ^. Using 
the explicit results of that reference, it is straightforward 
to ch eck t hat the model indeed obeys the Ward identity, 
Eq. (3.10).cZI This justifies, a posteriori, the approxima- 



tions that have entered our derivation. We will further 
discuss the merits and limitations of the model in the 
next section. 



V. DISCUSSION 

We finally discuss our results, and the procedures used 
to derive them. 



A. Role of the phase space decomposition 

Let us start wi th a discussion of the phase space decom- 
position in Sec. II A , which writes the interaction term 
(and also the disorder term) in the action as a sum of 
terms of different structures, with a cutoff to avoid double 
counting. As a result, the definition of the effective ac- 
tion contains this cutoff, which is a priori unspecified. For 

a (1,2,3) 



instance, without the cutoff A the three terms S'- 



Eqs. (|2.2|), would all be equal and equal to S^^. A super- 
ficial consideration might conclude that this decomposi- 
tion of the action introduces an unnecessary ambiguity 
into the theory. In fact, however, the phase space de- 
composition is necessary in order to derive a theory that 
allows for a well-behaved perturbation theory. This can 
be explained most easily by using the disorder term ^dis, 
Eqs. ( ^.11 ), as an example. As explained in I, the two 



(1 2) 

contributions ^^js result from a phase space decompo- 
sition analogous to the one performed on the interaction, 
and the matrix Q(x) is therefore to be understood as 
containing only Fourier components with wavenumbers 



|k| < A. Now suppose we had not performed the decom- 

(2) 

position. Then ^dis would consist of A]^^!. only, with Trei 
replaced by 2rrei, and Q(x.) containing all Fourier com- 
ponents. The saddle-point Green function for this action 
would then contain a disorder part of the self-energy that 
is half the Born value. In a perturbation expansion in 
powers of the disorder, higher orders would then have 
to make up for the missing factors of 2 at zeroth order, 
i.e. the perturbation expansion would be singular. This 
is precisely whatrJiappens in standard many-body per- 
turbation theoryJ3 where singular integrals make it im- 
possible to easily determine the order of a contribution 
from its diagrammatic structure, and infinite resumma- 
tions are in general necessary to obtain all contributions 
of a given order. In contrast, perturbation theory for 
the nonlinear a model is much better behaved, with the 
number of loops determining the order to which a given 
diagram contributes. This is a consequence of a judicious 
choice of the starting point for the loop expansion, which 
in turn depends crucially on the phase space decompo- 
sition. The fact that the theory depends on a cutoff is 
the price paid for the controlled nature of perturbation 
expansions. 

Similarly, the phase space decomposition performed on 
the interaction part of the action allows us to perturba- 
tively control the c ompl icated contribution W^^^ to the 



Ward identity, Eq. ( 3.11 ). This guarantees that the non- 
linear a model correctly reproduces perturbation theory 
in the interaction to first order, as is well kaowrufrom 
comparing results obtained by either method. Ea'tafl 

Conversely, the above discussion makes it clear that 
the nonlinear a model, owing to its derivation, is per- 
turbative in nature with respect to the electron-electron 
interaction. Indeed, its restriction in that respect is more 
serious than with respect to disorder: Since the loop ex- 
pansion is an expansion in powers of the disorder, going 
to higher order in perturbation theory will always include 
higher order disorder effects. With respect to the inter- 
action, the analogous statement is not true, since some 
effects of higher order in the interaction are left out of the 
model, although the loop expansion resums cetain classes 
of interaction terms to all orders. Of course, a complete 
renormalization of the model would in principle supply 
all of the effects that might have been left out of the bare 
model, but this will in general not be captured by the 
standard perturbative renormalization based on low or- 
ders in the loop expansion. This is a restriction that is 
important to keep in mind in the context of discussions 
about possible exotic effects of a strong effective interac- 
tion, like e.g. a metallic non- Fermi liquid ground state. It 
may also be relevant for understanding the observationtj 
that renormalization group calculations based on the a 
model in high dimensions (d > 6) reveal relevant terms of 
a structure that is not seen in low-order 2 + e expansions. 
Physically, the 'standard' generalized a model approach 
is valid if the physics one is interested in is determined by 
the two-particle diffusive modes. If, for example, there 
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were also soft single-particle excitations, then in general 
the coefficients in the a model would be singular and this 
approach would break down. 

Finally, we note that already in the Gaussian approxi- 
mation the perturbative nature of the a model approach 
with respect to interactions is apparent. In I we pointed 
out that before the massive modes were integrated out, 
the Gaussian field theory explicitly contained the Stoner 
theory for ferromagnetism. However, after the a model 
approximation was made, the interaction terms that lead 
to the Stoner theory were absent. 



B. Screening, and the disordered Fermi-liquid fixed 
point 

A characteristic feature of the long-ranged Coulomb in- 
teraction is that it leads to the system being incompress- 
ible: The wavenumber depended thermodynamic deriva- 
tive (i9n/9yLt)(k) = x(k, — 0), which is proportional to 
the compressibility, is for small wavenumbers given by 

TT^, ' (5-1) 



with K the screening wavenumber, and {dn/ d^)sc = 
Xsc(k — > 0,51 = 0) the screened density susceptibility, 
which is a nonzero number. This structure follows from, 
and is co ntrolle d by, the Ward identity, as can be seen 
from Sec. Ill B above. It is instructive to check explic- 



itly that the nonline ar a model respects the compress- 
ibility sum rule, Eq. (5.1). Within the framework of the 
a model, one hasB 



(5.2) 



with H as defined after Eq. (1.7), and ^^(k) = 
— 7rr(*^(k)/8. We see that the a model indeed respects 
the compressibility sum rule, with the Gaussian approx- 
imation for the screened compressibility, {dn/d^)sc ~ 
Np. Notice— that in the model originally proposed by 
Finkerstein,E3 Fermi liquid corrections had been put in 
to make {dn/dy)sc the exact screened density suscepti- 
bility. These are missing here, since we have integrated 
out the massive modes, which account for the screening, 
in a Gaussian approximation. We emphasize that, while 
it is of course always possible to put in Fermi liquid cor- 
rections by hand, nothing is really gained by doing so: 
Such a procedure just amounts to a partial resummation 
of some terms that are of higher order in the interac- 
tion, which docs not change the fact that the effective 
theory has a perturbative character with respect to the 
interaction, as was discussed in Sec. V A above. Further- 
more, the point of any effective theory is that it correctly 
captures the structure of the full theory, while the coeffi- 
cients can be represented by some approximation in the 
bare theory. Upon renormalizing the bare effective the- 
ory, the coefficients will get rcnormalized by fiuctuation 
effects. 



C. Renormalization group properties of the effective 
field theory 

We finally mention that t he r enormalization group 
properties of the a model, Eq. ( |4.7D , are well known. The 
theory possesses a critical fixed point that-dsscribes an 
Anderson-Mott metal-insulator transition.E2ll3 Also, due 
to the short-ra nged nature of the effective, screened, in- 
teraction (Eq. (4.6)), the discussion of the stable Fermi- 
liquid fixed point given in Sec. 111.B.2 of I still applies. 
The Fermi-liquid ground state is stable for c? > 2 in the 
presence of quenched disorder, and for d > 1 in the 
clean limit. The corrections to scaling yield the weak- 
localization nonanalyticities and their clean counterparts 
as discussed in I, modified by the log-squared singular- 
ities in the density of states and the sound attenuation 
that are induced by t he Coulomb interactionoQ as men- 
tioned after Eq. (3.10) above. We emphasize again, how- 
ever, that due to the perturbative nature of the effective 
field theory, our considerations do in no sense constitute 
a proof that the Fermi-liquid fixed point will be stable 
for arbitrary strengths of the bare interaction constant. 
What we have shown is that the fixed point is perturba- 
tively stable for weak interactions. It is easy to see that 
an interaction that is of longer range than the Coulombic 
one, ti(x) ^ l/|x|" with a < 1, will destroy the screen- 
ing process, and hence lead to a relevant operator that 
destroys the Fermi liquid fixed point, at least close to 
its lower critical dimension. If an effective interaction of 
such long range were generated by the renormalization 
group acting on the a model, or if it were present in a 
the bare action for a different effective theory that is not 
subject to the perturbative restriction of weak interac- 
tions, then this could lead to a non- Fermi liquid ground 
state. These points may be important in the context of 
the ongoing discussion about possible non-Fetpii liquid 
ground states in 2-d (clean) electron systems.Eil 
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APPENDIX A: PROPERTIES OF THE 
Q-MATRICES 

Here we list again the symmetry properties of the Q- 
matrices that were defined in I, and add some additional 
remarks. 
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Q is self-adjoint under the adjoint operation Q"*" 

C, with C^^ = Snm <5q/3 «n S2 , 



(Ala) 



In addition, the hermitian conjugate of Q is related 
to Q by,E3 



= (t3®so) rgr-^ (rg^so) 



(Alb) 



where the similarity transformation denoted by F has the 
property 



(rQr-i)„„. = 



-n — l.—7n — l 



(Ale) 



We now expand our matrix fields in the spin-quaternion 
basis defined after Eq. (|2.9|), 



Qi2(x) 

Al2(x) 



3 

E 

3 

E 

r,i=0 



(t^ «) Si) *Qi2(x) 

(r^ «) Si) ;A^2(x) 



(A2a) 
(A2b) 



where again 1 = (ni,ai), etc. In this basis, we have the 
following symmetry properties. 



rVl2 






?21 ' ('' = 0,3) , 


(A3a) 


rQl2 






;Q2i , (r = 0,3; ^ = 1,2,3) , 


(A3b) 


rQl2 


= ° 


921 


, (r = l,2) , 


(A3c) 


rQl2 




rQ21 


, (r = 1,2; z = 1,2,3) . 


(A3d) 



Toge ther with the behavior under hermitian conjugation, 
Eq. (Alb), this further implies 



rQl 



r^-ni-l,-n2-l 



(A4) 



for all i and r. Analogous relations hold for A by virtue 
of the linear coupling between Q and A. 



We further note that as a result of Eqs. (A3d) and 

^ 1,2.3, 



( A4) the matrix elements Qw ™ the particle-particle 
spin-triplet channel are real (not imaginary, as erro- 
neously stated in I). As a result, the Gaussian theory 
(Eqs. (2.36) in I) is formally unstable in that channel, 
and the formally diverging Gaussian integral needs to be 
interpreted. We have ascertained that a rotation of the 
relevant Q- integration contour onto the imaginary axis, 
which effectively makes 1'2'^Qj^^ imaginary, provides an 
interpretation that guarantees agreement with the, well- 
known results of conventional perturbation theory.H Also 
note that within the model, the matrix elements of Q ia. 
the particle-particle sp in-tr iplet channel are imaginary.cl 
This follows from Eq. ( Ala ) in conjunc tion w ith the her- 
miticity of Q within the a model, Eq. (4.4a). 
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